JLAB-THY-08-917 



Asymmetries involving dihadron fragmentation functions: from DIS to e+e 

annihilation 

Alessandro Bacchetta/'Q Federico Alberto Ceccopieri,^'[t| Asmita Mukherjee,^'[l and Marco Radici^'[§| 

'Theory Center, Jefferson Lab, 12000 Jefferson Ave, Newport News, VA 23606, USA 
^ Dipartimento di Fisica, Universita di Parma, 
and INFN, Gruppo Collegato di Parma, 1-43100 Parma, Italy 
''^Physics Department, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India 
^Istitwto Nazionale di Fisica Nucleare, Sezione di Pavia, 1-27100 Pavia, Italy 

Using a model calculation of dihadron fragmentation functions, we fit the spin asymmetry recently 
extracted by HERMES for the semi-inclusive pion pair production in deep-inelastic scattering on a 
transversely polarized proton target. By evolving the obtained dihadron fragmentation functions, 
we make predictions for the correlation of the angular distributions of two pion pairs produced in 
electron-positron annihilations at BELLE kinematics. Our study shows that the combination of 
two-hadron inclusive deep-inelastic scattering and electron-positron annihilation measurements can 
provide a valid alternative to Collins effect for the extraction of the quark transversity distribution 
in the nucleon. 
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I. INTRODUCTION 



Dihadron fragmentation functions (DiFF) describe the hadronization of a quark in two hadrons plus other unob- 
served fragments. In their simplest form, they represent the probability that at some hard scale a parton hadronizes 
in two hadrons with fractional energies zi and Z2- They were introduced for the first time when studying the 
e^e^ — !■ hih2X process in the context of jet calculus They are in fact necessary to guarantee the factorization of 
all collinear singularities for such a process at next-to- leading order (NLO) in the strong coupling constant 

In experiments, often not only the fractional energies of the two hadrons are measured, but also their invariant 
mass Mh (see, e.gj, Refs. [1, U (H). Hence, it is useful to introduce extended DiFF (in analogy with extended 
fracture functions [6|), which are explicitly depending on Mh- Their definition and properties have been analyzed 
up to subleading twist 0, Q . Their evolution equations are known Q and presently solved in the leading logarithm 
approximation (LL), and there are valid arguments to assume that they can be factorized and are universal, similarly 
to what happens for extended fracture functions In fact, (extended) DiFF can appear also in two-particle- inclusive 
deep-inelastic scattering (SIDIS) and in hadron-hadron collisions. 

DiFF can be used as analyzers of the polarization state of the fragmenting parton [HI ■ Because of this, they 
have been proposed as tools to investigate the spin structure of the nucleon, in particular to measure the transversity 
distribution hi of a parton q in the nucleon N (see Ref. [l3| for a review). The hf, together with the momentum 
fi and helicity gf distributions, fully characterizes the (leading-order) momentum/spin status of q inside N, if quark 
transverse momentum is integrated over. Transversity is a chiral-odd function and needs to appear in a cross section 
accompanied by another chiral-odd function. The simplest example is the fully transversely polarized Drell-Yan 
process, where /if appears multiplied by its antiquark partner [i4|. Although this process is theoretically very 
clean, it appears to be experimentally very challenging llSl , at least at present facilities (the same finding is confirmed 
for proton- proton collisions leading to prompt photon [l6| and semi- inclusive pion production [13 )• 

An alternative approach, so far the only fruitful one, is to turn to SIDIS and measure the correlation between the 
transverse polarization of the target and the transverse momentum of the final hadron, which involves a convolution 
of hi with the chiral-odd Collins fragmentation function H^'^ |T8| . The resulting asymmetry has already been 
measured at HERMES [H, [l^, and at COMPASS [U [H, [sl] . The knowledge of the Collins function is required 
to extract the transversity distribution. This can be obtained through the measurement of azimuthal asymmetries 
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in e+e" iT~^n~X with almost back-to-back pions [2^. The BELLE collaboration at KEK has measured this 
asymmetry [2^, making the first-ever extraction of hi possible from the global analysis of SIDIS and e~^e~ data [26| . 

At present, large uncertainties still affect this analysis and the resulting parametrization of hf. The most crucial issue 
is the treatment of evolution effects, since the BELLE and the HERMES/COMPASS measurements happened at two 

very different scales: ~ 100 and (Q^) = 2.5 GeV^, respectively. Both hi (8) and (g) Hf convolutions involve 



transverse-momentum dependent functions (TMD) [27l . |28| whose behaviour upon scale change should be described 
in the context of CoUins-Soper factorization |29l. [30| (see also Refs. [3l|, [s^). However, the global analysis of Ref. [2^ 
neglects any change of the partonic transverse momentum with the scale leading to a possible overestimation of 
/ti [28l. [33I. [3^ . It would be desirable, therefore, to have an independent way to extract transversity, involving collinear 
fragmentation functions. Here, we consider the semi-inclusive production of two hadrons inside the same current jet. 

As already explained, the fragmentation q (7r+7r^)X is described by an (extended) DiFF. When the quark is 
transversely polarized, , a correlation can exist between its transverse polarization vector and the normal to the plane 
containing the two pion momenta. This effect is encoded in the chiral-odd polarized DiFF via the dependence 
on the transverse component of the pion pair relative momentum Rt [3] . The function H ^'^ can be interpreted as 
arising from the interference of (tt+tt^) being in two states with different angular momenta |35l. [36l. [stI. [s^ . Since the 
transverse momentum of the hard parton is integrated out, the cross section can be studied in the context of collinear 
factorization and its polarized part contains the factorized product h\H^'^ jSTLISOj. The HERMES collaboration has 
recently measured such spin asymmetry using transversely polarized proton targets [i^; the COMPASS collaboration 



performed the same measurement on a deuteron target [41| and should soon release data using a proton target. 

Similarly to the Collins effect, the unknown Hf^ has to be extracted from electron-positron annihilation, specifically 
by measuriiig the angular correlation of planes containing two pion pairs in the e^e~ (7r"'"7r~)icti (7r''~7r~)jot2Ar 
process [l^. |42|. The BELLE collaboration is analyzing data for this angular correlation [H, |43|- Therefore, it seems 
timely to use available models for extended DiFF to make predictions for the e+e" azimuthal asymmetry at BELLE 
kinematics. Since evolution equations for extended DiFF are available at NLO @, at variance with the Collins effect 
the asymmetries with inclusive hadron pairs in SIDIS and e^e^ can be correctly connected when the scale is ranging 
over two orders or magnitude. Therefore, the option of using the semi-inclusive production of hadron pairs inside 
the same jet seems a theoretically clean way to extract transversity js^l- Finally, we point out that pair production 
in polarized hadron-hadron collisions allows in principle to "self-sufhciently" determine all the unknown DiFF and 
hi [4^ . The STAR collaboration has recently presented data on this kind of measurement [i^ . 

The paper is organized as follows. In Sec. |TT1 using the model calculation of DiFF from Ref. [Ill, we fit the 
spin asymmetry recently extracted by HERMES for the SIDIS production of (tt+tt^) pairs on transversely polarized 
protons ^4(7|. In Sec. IIIIl we describe how we calculate the evolution of the involved extended DiFF starting from 
the HERMES scale up to the BELLE scale. In Sec. IIVI we illustrate the predictions for the correlation of angular 
distributions of two pion pairs produced in e^e~ annihilations at BELLE kinematics. Finally, in SeclVlwe draw some 
conclusions. 

II. FIT TO DEEP INELASTIC SCATTERING DATA 

We consider the SIDIS process e{l) + N'^ {P) e{l') + 7r+(Pi) -I- tt^ {P2) + X, where P is the momentum of the 
nucleon target with mass M, I, I' are the lepton momenta before and after the scattering and q = I — I' is the space- 
like momentum transferred to the target. The final pions, with mass m,r = 0.14 GeV and momenta Pi and P2, have 
invariant mass Mh (which we consider as much smaller than the hard scale — —q^ > of the SIDIS process). We 
introduce the pair total momentum Ph = Pi + P2 and relative momentum R = (Pi — P2)/2. Using the traditional 
Sudakov representation of a 4- momentum a in terms of its light-cone components — {aP ± a?')/^/2 and transverse 
spatial components a^, we define the light-cone fractions x = jP^ and z = P^ /k~ , where p and k — p + q are the 
momenta of the parton before and after the hard vertex, respectively. 

In this process, the following asymmetry can be measured (for the precise definition we refer to Refs. [38l[40|). 



Ul \ Tfl l tlJ T I ,,.2 /o I ..2 -.2 /-I / - 2 \ 71 ,r ^2 -PQ / ^\ n l\/f2\ ' V / 



^^(1+7? (l + i^) \R\ E.elhUx) <7(.,M2) 



(1 + el Ilix) Di^,{z, Ml) 



where y = P ■ q/P ■ I is related to the fraction of beam energy transferred to the hadronic system, 7 — 2AIx/Q, ff 
and hf are the unpolarized and transversely polarized parton distributions, respectively, and 
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The spin asymmetry ([l} is related to an asymmetric modulation of pion pairs in the angles 4>g and which 
represent the azimuthal orientation with respect to the scattering plane of the target transverse polarization and of 
the plane containing the pion pair momenta, respectively (see Ref. [s^ for a precise definition, which is consistent 
with the Trento conventions [47|). 

The polar angle 6 describes the orientation of Pi, in the center-of-mass frame of the two pions, with respect to the 
direction of Ph in the lab frame. DiFF depend upon cos6' via the light-cone fraction ^ = 2R~ /P^ — 2cos9\R\/Mh, 
which describes how the total momentum of the pair is split between the two pions (38j . DiFF can be expanded 
in terms of Legendre polynomials of cos 9 and the expansion can be reasonably truncated to include only the s and 
p relative partial waves of the pion pair, since their invariant mass is small (typically Mh ^ 1 GeV). At leading 
twist and leading order in a^, the spin asymmetry of Eq. ((T|) contains only D^^ = Dl^ + D\ which includes the 

diagonal pure s- and p-wave contributions, and H^^^, which originates from the interference between them [ssl [39|. 
Subleading- twist terms have different azimuthal dependences [sj. 

The polarized DiFF H^^^ represents the chiral-odd partner to isolate the transversity distribution h\ in the spin 
asymmetry of Eq. ([T]). It describes the interference between the fragmentations of transversely polarized quarks 
into pairs of pions in relative s and p waves I39l . Together with it was analytically calculated in a spectator- 

model framework for the first time in Ref. [37| , and later in a refined version |38|. The analytical expressions 
for = -Di oo and = H^^^^ can be found in Eqs. (23) and (26) of Ref. [33, respectively. The model 

parameters were fixed by adjusting ^ to the output of the PYTHIA event generator tuned to the SIDIS kinematics 
at HERMES their values are listed in Eqs. (32-35) of Ref. 38]. Note that the spectator model by construction 
gives Di,, = D, , = D,^, = D^ ^ and H^:: = -H^' = -i/^? = ^1 ' 

The calculated spin asymmetry follows the same trend of the data. In particular, the shape of the invariant mass 
dependence is dominated by a resonance peak at Mh ~ rrip, which is due to the interference between a background 
production of pion pairs in s wave and the p— wave component dominated by the decay p — *■ tttt of the p resonance. 
Similarly, the model displays also another broad peak at ~ 0.5 GeV due to the lo — > tttttt resonant channel. Both 
predictions and data show no sign change in A!^^'^'^^'^^^ as a function of Mh around Mh ~ mp, contrary to what 
was predicted in Ref. [s^. However, the results of Ref. [Ill systematically overpredict the experimental data at least 
by a factor of two [i^ . 

To correctly reproduce the size of the asymmetry, we multiply the model prediction of H^'^^ by an extra parameter 
a, while we use the model prediction for ^ without further changes, since its parameters have been already fitted 
to reproduce the unpolarized cross section, as predicted by PHYTIA. We use also the GRV98 LO parametrization 
for fl at the HERMES scale = 2.5 GeV'^. For we use the recently extracted parametrization from Ref. ^48J, 
whose central value is basically the same as the former one from Ref. [IHl in the region x < 0.2 of interest here. The 
asymmetry is calculated by averaging the numerator and denominator of A^^^'^^''^'^^^ '^'"^ in each experimental bin, 
while integrating in turns all remaining variables in the ranges 

0.023 < a; < 0.4, < z < 0.99, 0.5 GeV < Ah < 1 GeV. (3) 

The variable y is always integrated in the range 

y,„i„ = Max[0.1, QlJ{x{s - M^)), {Wl,^ ~ M')/{{1 - x){s - M^))] , (4) 

with s = 56.2 GeV^ and W^^^ = 4 GeV^. 

The best value for the parameter a is found by means of a fit; the results are shown in Fig. [TJ We took into 
consideration the experimental errors by adding in quadrature the statistical and systematic errors (the error bands in 
Fig. [1] represent such a sum). We did not include the small theoretical errors coming from the uncertainty on the other 
model parameters of the fragmentation functions, from the uncertainties of the parton distribution functions, and 
from the choice of factorization scale. The best-fit value of our reduction parameter turns out to be a = 0.32 ± 0.06 
corresponding to x^/d.o.f. = 1.24. In summary, the model calculation of the H^^^ function has to be reduced by a 
factor 3 to reproduce the HERMES data, if the transversity from Ref. |4^ is used. 

III. EVOLUTION OF DIHADRON FRAGMENTATION FUNCTIONS 

In order to predict the azimuthal asymmetry in the distribution of two pion pairs produced in e^e~ annihilation, 
we need to evolve the DiFF D^ ^ and TJ^^^^^ from the HERMES scale to the BELLE scale. 



DiFF usually depend on z, C = 2cos9\R\/Mfi [or, alternatively, on zi = z{l -\- C)/2, Z2 = z{l — C)/2], and on i?; 
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FIG. 1: The spin asymmetry for the semi-inclusive production of a pion pair in deep-inelastic scattering on a transversely 
polarized proton, as a function of the invariant mass Mh of the pion pair, of the light-cone momentum fraction x of the initial 
parton, of the energy fraction z carried by the pion pair with respect to the fragmenting parton. Data from Ref. [i^. The 
uncertainty band is a fit to the data based on the DiFF spectator model of Ref. [38| and on the hi parametrization of Ref. [48l |. 



which is connected to the pair invariant mass by [9(] 

r Ml Ml Mi] 



2 (Pit - P2t) ziZ2 
rirp — ■ — 



Zl + Z2 



Zl + Z2 Zl Z2 



The further dependence on the scale of the process is described by usual DGLAP evolution equations; at LL, they 
read Q 



where P{u) are the usual leading-order splitting functions [i^ . A similar equation holds for H'^ involving the splitting 
functions 5P(u) for transversely polarized partons [13, HH (see also the Appendix of Ref. 0, for convenience). 

The same strategy can be applied to study evolution of single components of extended DiFF in the expansion in 
relative partial waves of the pion pair. In fact, Eq. ([6]) can be rewritten as 



d 



dlogQ 



D^{z, C, Ml Q^) = j" ^ (i, C, Ml Q^) P,. 



(") ■ (7) 



Note that the evolution kernel affects only the dependence on z, leaving untouched the dependence on C,. That is, it 
affects the dependence on the fractional momentum of the pion pair with respect to the hard fragmenting parton, but 
not the dependence on the nonperturbative processes that make the fractional momentum split inside the pair itself. 
The net effect is that extended DiFF display evolution equations very similar to the singlc-hadron fragmentation case. 
Using the above identity — 2 cos 0\R\/Mh, we can again expand both sides of Eq. ([7]) in terms of Legendre functions 
of cosO and apply the evolution kernel to each member of the expansion. By integrating in dcosO both sides we come 
to the final result 

-^D,Jz,MlQ') = ^4^ ['^D,^^,(^^,MlQ') P,,,{u), (8) 



dlogQ2 ^.9^ ' ' 2n J, u \u 



that involves the DGLAP evolution of the single diagonal component Di ^ Df + related to the pure s and 



p relative partial waves of the pion pair. Analogously, we can get an evolution equation similar to Eq. ([8]) for H^^^ 
provided that P{u) is replaced by SP{u). 

Equation ([5]) shows that also the dependence on the pair invariant mass M^ is not affected by the evolution kernel, 
as is reasonable, since Mh is a scale much lower than Q^. However, in order to get the Mh dependence at a different 
scale Q'^ 7^ it is important to completely integrate away the z dependence. Usually, experimental phase spaces are 
limited by the geometry of the apparatus and, in this case, the integration in dz is performed in the interval [^min, 1] 
with Zmin 0. In Fig. H we show D^^iMh) for the up quark at the HERMES scale = 2.5 GeV^ (dot-dashed 
line) and at the BELLE scale of = 100 GeV^ (solid line). In the left panel, results are obtained using Zy^in — 0.02, 
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FIG. 2: The unpolarized extended DiFF Dx.u(Mh,Q^) in arbitrary units, after integrating the z dependence away in the 
interval [0.02, 1] (left panel) and [0.2, 1] (right panel). Dot-dashed line for = 2.5 GeV^ at HERMES, solid line for = 100 
GeV^ at BELLE (see text). 

in the right panel with Zmin = 0.2 as in the BELLE setup. Since the DGLAP evolution shifts the strength at lower 
z for increasing Q^, cutting the z phase space from below makes the final result miss most of the strength and, 
consequently, display a reduced ^{Mh). The apparent (and contradictory) effect of perturbative evolution on the 
dependence upon the nonperturbative scale Mh in the right panel is actually spurious, and it disappears as soon as 
the phase space for z integration is properly enlarged to include the lowest z for Zmin — 0, as shown in the left panel. 
When extracting azimuthal (spin) asymmetries, it is, therefore, crucial to keep in mind these features to estimate the 
effect of evolution. 

As a last general comment, we stress that the analysis of evolution is facilitated by the fact that azimuthal asym- 
metries based on the mechanism of dihadron fragmentation can be studied using collinear factorization. This feature 
makes them a cleaner observable than the Collins effect, from the theoretical point of view. 

In the next section, we compute azimuthal asymmetries for two pion pairs production in e^e^ annihilations using 
evolved extended DiFF. The goal is to make model predictions at BELLE kinematics, and, at the same time, to 
estimate the evolution effects, both the pure one on the z dependence and the spurious one on the Mh dependence, 
due to the limited experimental phase space. 
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FIG. 3: Momenta and angles involved in the process e^e (tt+tt )jeti (7r'''7r )jet2X. 
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IV. PREDICTIONS FOR ELECTRON-POSITRON ANNIHILATION 

For the process e^e^ —f (7r+7r^)jcti(7r^7r^)jot2^, the momentum transfer q = I + I' is time-hke, i.e. ^ > 0, 
with 1,1', the momenta of the two annihilating fermions. We have now two pairs of pions, one originating from a 
fragmenting parton and the other one from the related antiparton. Therefore, we will use in analogy to Sec. |TT] the 
variables (p^^, 9, Pi, P2, Ph, R, Mh, z, (, for one pair, adding the variables 0^, 9, Pi, P2, Ph, R, M^, 'z, Q, for the 
other pair. Since we assume that the two pairs belong to two back-to-back jets, we have Ph ■ Ph ~ ■ The momenta 
and angles involved in the description of the process are depicted in Fig. [S] The azimuthal angles 0^ and 0^ are 
defined by 



X Ph) ■ Rt ( le+ X Ph Rt X Ph 

X Ph) ■ Rt\ xPh\' \Rtx P, 



'h 



- [h+ X Ph) ■ Rt f{le+ X Ph) {Rt x Ph) 

f>R — = — arccos -r, 7:rT • ^= 

'^'^ X Ph) ■ Rt\ \\h+xPh\ \RtxP, 



(9) 



h\ 



where /g+ is the momentum of the positron, and Rt, Rt indicate the transverse component of R, R with respect 
to Ph, Ph, respectively. They are measured in the plane identified by 1^+ x Ph and {le+ x Ph) x z, with z || —Ph 
in analogy to the Trento conventions [13]; this plane is perpendicular to the lepton plane identified by ie+ and Ph 
(see Fig.[3l). Note that the difference between as defined in Eq. ([9]), and the azimuthal angle of Rt, as measured 
around Ph, is a higher-twist effect. The invariant y — Ph - IjPh ■ q is now related, in the lepton centcr-of-mass frame, 
to the angle 92 = arccos(Z^+ • Ph/{\l^+\ \Ph\)) by y = (1 + cos6'2)/2. 

Starting from Eq. (21) of Ref. [4^, we define the so-called Artru-CoUins azimuthal asymme- 

try A(cos ^2, -2:, z, M^, M^^) as the ratio between weighted leading-twist cross sections for the e^e" 

2 

(7r+7r^)joti(7r"'"7r~)jet2-'^' process, once integrated upon all variables but cos 92, z,z, M'^, M^. We define the weighted 
cross section as 

(w) = f dCd(^ I d4>j^d(t)fi I dqT w _ — ^^2 = • (10) 

J J J d cos 92 dzdzdCdCdMldMf^d(j)j^d(t)j^ dqT 

The weight of the numerator in the asymmetry is cos((/)^ + (t>fi), for the denominator it is just 1. Elaborating upon 
the work of Ref. [13] ^, we recall the change of variables ( = 2 cos 6'|il|/Af;i (and similarly for () and we perform an 
expansion in terms of Legendre functions of cos 9 (and cos 9) by keeping only the s- and p-wave components of the 
relative partial waves of the pion pair. By further integrating upon dcos9 and dcos9, we deduce the analogue of 
Eq. (32) of Ref. [i^] for the specific contribution of s and p partial waves to the Artru-CoUins azimuthal asymmetry, 
namely 

Aicos92, z,Ml^, Ml) ^^^'^^-^^-^^ 



1 

sin^ 92 7:^ \R\ \R\ H^;/{z, Ml) H-;(z, Ml 



1 + cos2 02 32 M,M, Y.,^lD,,^{z,Ml)D,^^{z,Ml) 

The extended DiFF ^ and -ff^^^ are the same universal functions appearing in the SIDIS spin asymmetry of 
Eq. ([1]). Hence, tuning model predictions for them at BELLE kinematics would help in reducing the uncertainty in 
the extraction of the transversity hi at the HERMES scale. In this strategy, a crucial role is played by evolution. 
At variance with the Collins effect, the dihadron fragmentation mechanism is fully coUinear, since only Rt matters 
and Ph± can be integrated. Hence, evolution equations for extended DiFF are easily under control, presently at LL 
level 0, and are represented by Eq. ^ and its analogue for i?^^^. 

In Fig.m the azimuthal asymmetry of Eq. (fTT|) is displayed as a function of Alh for the z bins [0.01,0.1], [0.1,0.2], 
[0.2,0.3], [0.3,0.4], [0.4,0.55]', and [0.55,0.75], after integrating upon the other variables 0.4 < M/^ < 1.2 GeV, 



^ With respect to Ref. [42l | , we use the modified definition ^ = 2^ — 1 , the integration measure of the parton transverse momentum reads 
z-^dkx, and the scaUng factor l/(Mi -|- M2) of the chiral-odd projections of the parton-parton correlator is replaced by l/M^ (and 
similarly for the antiparton correlator), in agreement with Refs. [Sa. |39(I . 
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FIG. 4: The azimuthal asymmetry for two pion pairs production in e e~ annihilation as a function of the invariant mass Mh 
of one pair for the indicated bins in its momentum fraction z. Notations as in Fig. (2] The uncertainty band around the sohd 
line originates from the fit error of Fig. [T] through error propagation. For each panel, the lower plot shows the modification 
factor of the final result because of DGLAP evolution (see text). 



0.2 < z < 0.9, and —0.6 < cos 6*2 < 0.9, according to the BELLE experimental phase space. In particular, according 
to Ref. [1^ for each bin we have assumed the following coefhcient 



wO.7. (12) 
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In each panel, the upper plot shows the A{zhin, M'^) at the HERMES scale = 2.5 GeV^ (dot-dashed line) and 
at the BELLE scale — 100 GcV^ (solid line), the latter being supplemented by the uncertainty band propagated 
from the SIDIS fit error of Fig. [TJ The lower plot shows 



A = 



A(zbin,M^,( 



100)-A(zbin,M2,Q2 ^2.5) 



A(zbin,M2,Q2 = 100) 



(13) 



i.e. the modification factor of the final result due to the evolution starting at the HERMES scale. 

Some comments are in order about Fig. [D First of all, the absolute size of the Artru-CoUins asymmetry is small, 
reaching at most the percent level (see Fig. [5]). However, it should be within the reach of the BELLE experimental 
capabilities, if compared with the corresponding Collins effect for two separated single-hadron fragmentations [25j . 
The error band originates from the uncertainty in the size of DiFF due to the fit of the SIDIS data for the spin 
asymmetry of Eq. ([1]). This error band is always much larger than the effects due to evolution. 

After the comments about Fig. [2l one would be tempted to attribute this sensitivity of the Mh dependence to the 
hard scale as coming from a spurious effect; indeed, in each panel of Fig. |4] the asymmetry is integrated in the 
indicated z bin, which is obviously just a small fraction of the available phase space. However, the asymmetry of 
Eq. pT|) is the ratio of two objects that behave very differently under evolution because of their kernels P{u) and 
SP{u), respectively. Hence, there is no fundamental reason to expect the pure Mh dependence of the asymmetry be 
preserved by DGLAP evolution, even after integrating upon the whole phase space of the other variables. Moreover, 
the moderate sensitivity of ^(zbin, Af^) to the hard scale of the process is the result of a compensation of two very 
large sensitivities both in the numerator and in the denominator, as is clear from Fig. [5l 
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FIG. 5: Numerator (left panels) and denominator (right panels) of the azimuthal asymmetry in the same conditions and with 
the same notations as in Fig. |4l for the indicated boundary bins in z. The result at = 100 GeV^ (solid line in the upper 
plot of each panel) is emphasized by the factor 10. 
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In Fig. [5l the numerator (left panels) and denominator (right panels) of the asymmetry pT|) are shown with the 
same notations as in Fig. 0] for the 0.01 < z < 0.1 (upper panels) and 0.55 < z < 0.75 bins (lower panels). The solid 
line, corresponding to the result at = 100 GeV^, is amplified by a factor 10. Therefore, the effect of DGLAP 
evolution is enormous, both in the numerator and in the denominator, where, in particular, it can reach a reduction 
factor of more than two orders of magnitude. Also the shape of the dependence is altered, making the more or 
less stable trend of A{zhin, M^) at different just a fortuitous case. 

In summary, even if DGLAP evolution of extended DiFF seems to mildly affect the predictions for the azimuthal 
asymmetry at BELLE, this small sensitivity arises from a dramatic compensation between big modifications in the 
numerator and in the denominator of the asymmetry. Therefore, it is wise to carefully consider such effect, because 
it could provide more sizeable modifications in other portions of the phase space. 

For sake of completeness, in Fig. [H] the azimuthal asymmetry pT|) is displayed as a function of z for the z bins 
[0.01,0.1], [0.1,0.2], [0.2,0.3] J0.3, 0.4], [0.4, 0.55], and [0.55, 0.75], with the same notations as in Fig. Hand again after 
integrating upon 0.4 < Mh,Mh < 1-2 GeV, and —0.6 < cos < 0.9. It shows a rising trend for increasing both z 
and z. The effect of DGLAP evolution is small and within 10%, except for the lowest z values. 

V. CONCLUSIONS 

In this paper, using the model calculation of extended dihadron fragmentation functions (DiFF) from Ref. (s^l 
we fitted the spin asymmetry recently extracted by the HERMES collaboration for the semi-inclusive deep-inelastic 
scattering (SIDIS) production of (tt+tt") pairs on transversely polarized protons [40^. Then, using the results of 
Ref. [Q] we calculated the evolution of extended DiFF at leading logarithm level, starting from the HERMES scale up 
to the scale of the process e+e" (7r+7r^)jcti (7r+7r~)jet2-''^ at BELLE kinematics. Finally, we made our predictions 
for the so-called Artru-CoUins asymmetry, describing the correlation of angular distributions of the involved two pion 
pairs. The BELLE collaboration plans to measure this asymmetry in the near future. 

The absolute size of the Artru-CoUins asymmetry turns out to be small, but it should be within reach of the 
BELLE experimental capabilities, if compared with the corresponding Collins effect for two separated single-hadron 
fragmentations [1^. The theoretical error band, originating from the uncertainty in the fit of the SIDIS data, is always 
larger than the effects produced by the evolution of DiFF. The latter seems to mildly affect the predictions for the 
azimuthal asymmetry at BELLE. Nevertheless, this small sensitivity arises from a dramatic compensation between big 
modifications in the numerator and in the denominator of the asymmetry. Therefore, it is wise to carefully consider 
such effect, because it could provide more sizeable modifications in other portions of the phase space. 

We stress that azimuthal asymmetries based on the mechanism of dihadron fragmentation can be studied using 
collinear factorization, which facilitates the analysis of, e.g., evolution. From the theoretical point of view, this feature 
makes them a cleaner observable than the Collins effect in single-hadron fragmentation, where transverse-momentum 
dependent (TMD) functions are involved, whose evolution is yet not taken into account. All this procedure would not 
be plagued by theoretical uncertainties about factorization and evolution of TMD parton densities, which currently 
affect the analysis of single-hadron fragmentation. As a consequence, the option of using the semi-inclusive production 
of hadron pairs inside the same jet seems theoretically the cleanest way to extract the transversity distribution hi , at 
present [34]. 

When BELLE data will be available, it will be possible to constrain extended DiFF on e+e^ data, to evolve them 
back to the HERMES scale, and to use them in the formula for the SIDIS spin asymmetry to directly extract hi. 
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